Quantum Critical Point at Finite Doping in the 2D Hubbard Model: A Dynamical 

Cluster Quantum Monte Carlo Study 
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We explore the Matsubara quasiparticle fraction and the pseudogap of the two-dimensional Hub- 
bard model with the dynamical cluster quantum Monte Carlo method. The character of the quasi- 
particle fraction changes from non-Fermi liquid, to marginal Fermi liquid to Fermi liquid as a 
function of doping, indicating the presence of a quantum critical point separating non-Fermi liquid 
from Fermi liquid character. Marginal Fermi liquid character is found at low temperatures at a 
very narrow range of doping where the single-particle density of states is also symmetric. At higher 
doping the character of the quasiparticle fraction is seen to cross over from Fermi Liquid to Marginal 
Fermi liquid as the temperature increases. 
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Introduction- The unusual properties of the hole 
doped cuprate phase diagram, including a pseudogap 
(PG) at low doping, and unusual metallic behavior at 
higher doping have lead many investigators to propose 
that there is a quantum critical point in the cuprate 
phase diagram at optimal doping. Different investiga- 
tors argue that the PG is related with the establishment 
of order P, 0, [3, 0, [Hli where the optimal doping is in 
the proximity of the quantum critical point (QCP) asso- 
ciated with this order [3| • Other investigators argue that 
the QCP is located at the transition from a non-Fermi 
liquid (NFL) to Fermi liquid (PL) ground state without 
the establishment of order in PG regionQ. 

In previous work employing cluster extensions of the 
Dynamical Mean Field, a PG region was found at low 
doping in the two-dimensional Hubbard model. It is 
characterized by a narrow gap-like feature in the single- 
particle density of states, a suppression of the spin sus- 
ceptibility, and a self energy of NFL character 
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In this paper, we employ significantly larger clusters 
than most previous studies, which affords us greater mo- 
mentum resolution. We investigate the single particle 
properties of the two-dimensional Hubbard model with 
the dynamical cluster approximation (DCA) [3, [l^ ■ We 
find further evidence for a QCP and are able to determine 
its character as the terminus of a Marginal Fermi Liquid 
(MFL) region, which separates a NFL PG region at low 
doping from a FL region at high doping. We present 
a comparative discussion of a few existing scenarios for 
quantum criticality in the context of our results. 

Formalism- We consider a 2D Hubbard Hamiltonian 



H 



-h.C 



(1) 



where t is the hopping matrix, cjj^(cio-) is the creation 
(annihilation) operator for electrons on site i with spin 
(T, and U is the on-site Coulomb repulsion which is taken 
to be three quarters of the bandwidth 8t). The 

hopping t is restricted to nearest neighbors (ij). 

We employ the DCA with a quantum Monte Carlo 
algorithm as the cluster solver. The DCA is a cluster 
mean-field theory which maps the original lattice model 
onto a periodic cluster of size Nc = L^. embedded in a 
self-consistent host. Spatial correlations up to a range 
Lc are treated explicitly, while those at longer length 
scales are described at the mean-field level. However 
the correlations in time, essential for quantum criticality, 
arc treated explicitly for all cluster sizes. To solve the 
cluster problem we use the Hirsch-Fye quantum Monte 
Carlo method 

and employ the maximum entropy 
method [l7| to calculate the real frequency spectra. 

A number of the normal state anomalies in the 
cuprates are describable by a MFL in proximity to a 
quantum critical point at finite doping^. The imagi- 
nary part of the MFL self energy has the form 



-a max (|a;| , T) 



(2) 



In contrast the imaginary part of the FL self energy has 
the form 



Y.'^l{uj) = -amax(w2^T2) 



(3) 



which would be expected to be valid at large doping and 
low temperatures. In the doping region beyond but near 
the QCP, the single particle properties of the model are 
observed to cross over from FL to MFL as the tempera- 
ture crosses Tx and the frequency w^, we find that this 
may be fit with the form 



-aujx max (|w| , T) for \uj\ > ujx or T > Tx 
-a max (cj^, T'^) for |cij| < oj^ and T < Tx 

(4) 
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Causality requires that a > 0, and the integrals over 
these forms for the self energy are cutoff at Wc which is 
of the order of the bandwidth W{= 8t). 

To compare our Matsubara frequency results to these 
forms of the real- frequency self energy [isj, we transform 
each of these forms to Matsubara frequency using the 
transform of the non-Hartree part of the self energy 



and then evaluate 



Zo(k) = (l-ImI](k,iwo)/'^o) 



(5) 



(6) 



where ldq = ttT is the lowest Fermion Matsubara fre- 
quency. For a well behaved self energy, limy^o ^o(k) = 
Z(k) is the quasiparticlc rcnormalization factor. For ex- 
ample, for the MFL, we find[li| 



Wo 



In 



■ tan 



1 



(7) 



While for the FL. wc find 



= I h 0.066235 — TT tan 

too TT \ T ttT / 

(8) 

when T < LOc- The crossover form is more complicated, 
but can be constructed from the same integrals used to 
derive Eqs. [7] and [5] 



ImI](iwo) 
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The parameters, a, lOc, Tx, and tOx are determined from 
a fit to the quantum Monte Carlo data. 

Results- We will present results for the model with 
t/ = 1.5 with bare bandwidth W = 2 setting the energy 
unit, and a 4 x 4 cluster. With this choice of U/t and 
cluster we are able to access low temperatures T ^ 0.01 
before the average sign of the sampling weight falls below 
0.05. The low energy scales in the problem are the an- 
tiferromagnetic exchange energy J near half filling, the 
PC temperature T* in the PC region, and the effec- 
tive Fermi energy Tx at higher doping. As described 
previously we can extract the effective near neighbor 
spin exchange energy Jej f from the spin excitation spec- 
trum at K = (0, tt). At this wavevector the magnon peak 
is expected at frequency oj ~ 2Jeff- From this analysis, 
we find that Jeff ~ 0.11 for N = 0.95 and iV = 1. The 
energy scales T* and Tx are extracted from fits to the 
data presented below. It is important to note that in 
each of these fits, we include data for T <C Jeff- 
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FIG. 1: Matsubara quasiparticle fraction Zo(k) versus tem- 
perature T evaluated with k on the Fermi surface along the 
(0, 1) direction for different fillings A'^ when U = 1.5 and the 
bandwidth W = 2. The lines represent fits in the region 
T < 0.3 to either the MFL form, Eq. for iV > 0.85 or the 
crossover form (X), Eq. [O] for A'' < 0.85. The arrows indicate 
the values Tx extracted from the crossover fits or T* (c.f. 
Figs. [3] and Note that the data for = 0.85 fits the MFL 
nearly perfectly, while the data for A'^ > 0.85 is poorly fit by 
the MFL for y <; y* which is too slow in temperature to pro- 
vide a good fit due to the formation of the pseudogap. The 
data for A'^ = 0.75 was also fit by the FL form, Eq.|8] however, 
the fit is clearly worse than that obtained by the crossover 
form. Inset: The ratio, Zqh/Zooi, is plotted as a function of 
temperature for different fillings. The ratio is essentially the 
same for all fillings at the QCP, indicating that Z is essentially 
isotropic, and becomes progressively more anisotropic as we 
dope into the PG region. 



The Matsubara quasiparticlc fraction is calculated 
with k on the Fermi surface defined by the maximum 
|Vn(k)| along the (1,1) and (0,1) directions. We will 
present detailed results and analysis for the latter only 
as we are interested in the crossover from PG to FL be- 
havior, and the PG is stronger along the (0, 1) direction. 
Zooi is shown in the main panel in Fig. [T] for U = 1.5 in 
units where W ~ 2 for different fillings. 

The low temperature Matsubara quasiparticle data 
changes character as the filling increases through 
N = 0.85. The data for A^ > 0.85 has negative cur- 
vature at all T. Whereas the data for < 0.85 has 
negative curvature at high T, a region of weak positive 
curvature is found at lower T. The change in curva- 
ture of the low temperature data for < 0.85 is easily 
understood as a crossover to a FL region. In a FL at 
zero temperature Zo_fl(0) = 1/(1 + 2auJc/T^) while for 
low T, Zofl{T) « Zofl(O) (1 + 3.099aZoFLiO)T). Since 
a^oFL(O) > 0, Zopl{T) at low T has a finite intercept 
and a linear region with positive initial slope indicative 
of FL formation. The next correction, of order T^, is 
small for Zqfl{0) « 0.6 so that Zofl{T) is a nearly lin- 
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FIG. 2; The single-particle density of states (DOS) when 
U = 1.5, W = 2, and T = 0.01429. Note that the low-energy 
DOS is roughly symmetric around lj — for the critical filling 
A'^ = 0.85 where Zo(k) fits a MFL form. This is consistent 
with the observation of rough particle-hole symmetry in the 
cuprates in the proximity of optimal doping [6[. 
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FIG. 3: The single-particle density of states in the pseudogap 
region for various temperatures with A*' = 0.95, U = 1.5, 
W = 2. Inset: The bulk, Q = 0, cluster susceptibility for the 
same parameters. The PG in the DOS begins to develop at 
roughly the same temperature T* which identifies the peak 
susceptibility. 



ear function when fit to the low T data. On the other 
hand, the MFL always has negative curvature, as can be 
seen from an expansion of Eqs. [7] and [5] to second order 
in T. So at the transition between FL and MFL, a region 
of positive curvature is found at T « Tx- The data for 
N < 0.85 is well fit by the crossover form posed above, 
but is poorly fit by a FL form over the fitting region (see 
e.g., the solid line fit to the N = 0.75 data). When the 
data at = 0.85 is fit with the crossover form for the 
Zqx{T), the fitting routine returns Ux = Tx = (within 
the precision of the fit and data), consistent with the for- 
mation of a MFL. So the solid line shown in the plot is 
a MFL fit. The MFL fit to the N = 0.85 data is very 
good. In fact the quality of this fit was better than that 
obtained for any of the fitting forms to any of the other 
data sets, despite the fact that the MFL form only has 
two adjustable parameters. In order to show that the 
conclusions from the above analysis are not specific to 
the direction (0,1), we plot the ratio, Zqh/^^ooi: in the 
inset of Fig. [T] as a function of temperature for different 
fillings. The ratio is seen to be essentially the same for 
all fillings at the QCP, indicating that Z is essentially 
isotropic at the QCP, and becomes progressively more 
anisotropic as we dope into the PG region. 

The data with = 0.85 in some other ways is special. 
For example, at this filling the low temperature single- 
particle density of states (DOS), which is plotted in Fig. [2] 
for several fillings, is peaked at zero frequency. At low en- 
ergies \ui\ ^ Jeff the A^ = 0.85 DOS is nearly symmetric 
around this point. This is consistent with the observa- 
tion of particle-hole (p-h) symmetry in the transport of 
the cuprates at optimal doping^. 

In order to characterize the region A^ > 0.85, the PG 



region, we also explored the temperature dependence of 
the DOS and the bulk (Q = 0) spin susceptibility of the 
cluster, as shown in Fig. [3] and its inset, respectively. We 
find a concomitant depression of the low energy DOS at 
temperatures below the peak in the susceptibility. The 
suppression of the susceptibility indicates the suppression 
of low energy spin excitations. Zqoi{T) in this region is 
well fit for T > T* (see Fig. IH) by the MFL form, but 
is poorly fit by the MFL form for T < T* , as shown in 
Fig.[TJ The MFL form changes too slowly with decreasing 
T, due to the formation of the PG for T <T*. 

The relevant temperatures near the QCP, Tx and T*, 
are shown in Fig. 2] The PG temperature was deter- 
mined from the peak in the susceptibility and the initial 
appearance of the PG in the DOS as shown in Fig. [3l 
and Tx from the fit to Eq. M 

We also explored the effect of larger clusters and U/t; 
however, these results are restricted by computational 
limitations including, especially, the minus sign problem. 
Despite this, some calculations were possible for Nc = 24 
site clusters. For the same parameters used above, we 
find that the PG temperature T* increases slightly near 
half filling. The critical filling where the Zqoi {T) is best 
fit by a MFL, the DOS is p-h symmetric at low frequen- 
cies remains between A^ w 0.85 to A/ « 0.86. Due to the 
minus sign problem, we did not have enough low temper- 
ature data to determine the crossover temperature pre- 
cisely for clusters with Nc = 24. We also explored larger 
values oi U = W = 2 for the Nc = 16 cluster (this data 
was produced in an earlier study). Again conclusions 
are limited by the minus sign problem; however, we find 
that the PG vanishes at roughly A^ = 0.78, consistent 
with optimal doping in the cuprates 0, and at this filling 
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FIG. 4: The pseudogap temperature T* , identified from the 
peak in the susceptibility and the emergence of the PG in the 
DOS shown in Fig. [3] and the FL to MFL crossover tempera- 
ture identified by fitting Eq.[9]to the Matsubara quasiparticle 
data shown in Fig. [T] 



the DOS is again p-h symmetric at low energies. 

Discussion- Several different scenarios have been pro- 
posed to explain the transition from PG to FL behavior 
and the strange behavior seen near optimal doping in the 
cuprates such as the MFL phenomena. 

An extension of the bond-solid scenario 0, 0, [l^ is 
consistent with our results. Here, a bond-solid is conjec- 
tured to coexist with the pseudogap in the underdoped 
region, and the QCP marks the doping where the tran- 
sition temperature to bond order vanishes. We suggest 
that the NFL behavior could be due to the scattering of 
quasiparticles from bond excitations. For dopings greater 
than, but near, the PG region, there would be remnant 
bond excitations. The low energy scale of these excita- 
tions will be cut off by the finite correlation length of 
the bond order, yielding a gap to bond excitations. So 
low energy quasiparticles may form a FL, while higher 
energy ones do not due to scattering from these bond ex- 
citations. In this scenario, the gap to bond excitations is 
proportional to Tx, which presumably will grow as the 
bond correlation length falls as the system is doped away 
from the QCP. One problem with this scenario is that it 
requires gapless bond excitations in the PG region, but 
due to the finite size of the cluster there should be a gap 
which scales with the cluster size. However, this gap may 
be small so its effects might not be seen at the tempera- 
tures we can access. 

The spectral- weight transfer scenario 0, [2lj also pro- 
vides a consistent interpretation. In a hole-doped Mott 
insulator, with doping a; = 1 — and large U /t, each 
doped hole yields two states immediately above the chem- 
ical potential. One comes from the lower and the other 
from the upper Hubbard band. When the system is 
doped so that the mmrber of low-energy states (i.e., not 



including those in the upper band) above and below 
the chemical potential are equal, so that 2x = 1 — x 
or X = 1/3, then p-h symmetry is obtained. For finite 
C//i, the critical doping for p-h symmetry is smaller 
This scenario also provides a general mechanism for the 
breakdown of the FL in the PG region [2^. In a FL 
there is a one-to-one correspondence between the par- 
ticles and quasiparticles which means that the number 
of ways of adding a particle at low energies equals the 
number of ways that one can add an electron to the un- 
occupied states. This correspondence breaks down in the 
PG region since the number of ways to add an electron 
is 2a;, while the number of ways one can add a particle 
is larger 0. Since such an argument relies on the asym- 
metry of adding a hole or an electron, it is clearly invalid 
when p-h symmetry is achieved for the low energy states. 
Thus this scenario explains the p-h symmetry of the DOS 
at the QCP, as found in Fig. H 

Our results differ from previous extended DCA results 
for the t-J model on a A'c = 4 cluster where a FL-FL 
crossover, with maximum scattering is found at a critical 
point not associated with MFL behavior It is not 
clear whether the model or the method is responsible for 
these differences, although note that the spectral- weight 
transfer arguments discussed above suggest FL physics 
for the t-J model even at small doping [2 2|. 

Conclusion- We investigate the Matsubara quasipar- 
ticle fraction on the Fermi surface and the PG of the 
two-dimensional Hubbard model. As a function of dop- 
ing, Zqqi{T) changes character. For doping beyond the 
critical point, as the temperature is lowered the curva- 
ture of Zooi(r) changes from negative to positive. This 
can be understood as a change from a MFL to a FL as T 
falls. At lower doping, the curvature is negative for all T, 
including T ^ Jeffi consistent with a NFL state. A PG 
is also found in the DOS and the bulk spin susceptibility 
at lower doping. At the QCP which separates these two 
regions, we find a MFL which is also found for T > T* 
and T>Tx. 
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